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£^ | In 1970s, Gutman introduced the concept of the energy S(G) for a simple graph 

^ \ G, which is defined as the sum of the absolute values of the eigenvalues of G. This 

■ graph invariant has attracted much attention, and many lower and upper bounds 

have been established for some classes of graphs among which bipartite graphs are of 
, particular interest. But there are only a few graphs attaining the equalities of those 

^) \ bounds. We however obtain an exact estimate of the energy for almost all graphs 

f-n | by Wigner's semi-circle law, which generalizes a result of Nikiforov. We further in- 

vestigate the energy of random multipartite graphs by considering a generalization 
of Wigner matrix, and obtain some estimates of the energy for random multipartite 
graphs. 
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o\ : 1 Introduction 
O 

Throughout this paper, G denotes a simple graph of order n. The eigenvalues Ai, . . . , A n of 
the adjacency matrix A(G) = (aij) nxn of G are said to be the eigenvalues of the graph G. 
In chemistry, the eigenvalues of a molecular graph has a closed relation to the molecular 
orbital energy levels of 7r-electrons in conjugated hydrocarbons. For the Hiichkel molecular 
orbital approximation, the total 7r-electron energy in conjugated hydrocarbons is given 
by the sum of absolute values of the eigenvalues of the corresponding molecular graph in 
which the maximum degree is not more than 4 in general. In 1970s, Gutman [U] extended 
the concept of energy $(G) to all simple graphs G, and denned that 

n 

*(G) = X>I. 
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where Ai, . . . , A n are the eigenvalues of G. Evidently, one can immediately get the energy 
of a graph by computing the eigenvalues of the graph. It is rather hard, however, to 
compute the eigenvalues for a large matrix, even for a large symmetric (0,l)-matrix like 
A(G). So many researchers established a lot of lower and upper bounds to estimate the 
invariant for some classes of graphs among which the bipartite graphs are of particular 
interest. For further details, we refer the readers to the comprehensive survey [TU]. But 
there is a common flaw for those inequalities that only a few graphs attain the equalities 
of those bounds. Thus we can hardly see the major behavior of the invariant S{G) for 
most graphs with respect to other graph parameters (|V((j)|, for instance). In this paper, 
however, we shall present an exact estimate of the energy for almost all graphs by Wigner's 
semi-circle law. Moreover, we investigate the energy of random multipartite graphs by 
employing the results on the spectral distribution of band matrix which is a generalization 
of Wigner matrix. 

The structure of our article is as follows. In the next section, we shall consider the 
random graphs constructed from the classical Erdos-Renyi model. The second model is 
concerned with random multipartite graphs which will be defined and explored in the last 
section. 

2 The energy of G n (p) 

In this section, we shall formulate an exact estimate of the energy for almost all graphs 
by Wigner's semi-circle law. 

We start by recalling the Erdos-Renyi model Q n (jp) (see [I]), which consists of all 
graphs with vertex set [n] = {1,2, ... ,n} in which the edges are chosen independently 
with probability p = p{n). Apparently, the adjacency matrix A(G n (p)) of the random 
graph G n (p) G G n (p) is a random matrix, and thus one can readily evaluate the energy of 
G n (p) once the spectral distribution of the random matrix A(G n (p)) is known. 

In fact, the study on the spectral distributions of random matrices is rather abundant 
and active, which can be traced back to [15]. We refer the readers to [21 El [7] for an 
overview and some spectacular progress in this field. One important achievement in that 
field is Wigner's semi-circle law which characterizes the limiting spectral distribution of 
the empirical spectral distribution of eigenvalues for a sort of random matrix. 

In order to characterize the statistical properties of the wave functions of quantum 
mechanical systems, Wigner in 1950s investigated the spectral distribution for a sort of 
random matrix, so-called Wigner matrix, 

X n := (xij), 1 < i,j < n, 

which satisfies the following properties: 

independent random variables with 

• the Xa's have the same distribution F±, while the x^s (i ^ j) have the same 
distribution F 2 ; 

• Vax(xy) = a\ < oo for all 1 < % < j < n. 
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We denote the eigenvalues of X n by Ai jn , A 2i n, • • • , X n>n , and their empirical spectral dis- 
tribution (ESD) by 

$x n (x) = - ■ #{A i>n | Xi n <x, i = 1, 2, . . . , n}. 
n 

Wigner [13], [H] considered the limiting spectral distribution (LSD) of X n , and obtained 
his semi-circle law. 

Theorem 2.1. Let X n be a Wigner matrix. Then 

lim $ n -i /2x (x) = a.s. 

n— »oo 

z.e. ; wrf/i probability 1, the ESD $ n -i/2 X „( a; ) converges weakly to a distribution as n 
tends to infinity, where has the density 



0( x ) = o 2y Aa l~ x2 1 kl<2<7 2 - 

Remark 2.1. One of classical methods to prove the theorem above is the moment 
approach. Employing the method, we can get more information about the LSD of Wigner 
matrix. Set fa = J xdFi (i = 1, 2) and X n = X n — fal n — ^(^n ~ In), where l n is the 
unit matrix of order n and J n is the matrix of order n in which all entries equal 1. It is 
easily seen that the random matrix X n is a Wigner matrix as well. By means of Theorem 
12.11 we have 

lim $ n -i/2x n (x) = $(x) a.s. (1) 

n^oo 

Evidently, each entry of X n has mean 0. Furthermore, one can show, using moment 
approach, that for each positive integer k, 

lim / i fc # n -i/2x (x) = / x k d§(x) a.s. (2) 

rwoc J J 

It is interesting that the existence of the second moment of the off-diagonal entries is 
the necessary and sufficient condition for the semi-circle law, but there is no moment 
requirement on the diagonal elements. For further comments on the moment approach 
and Wigner's semi-circle law, we refer the readers to the extraordinary survey by Bai [2] . 

We shall say that almost every (a.e.) graph in Q n {p) has a certain property Q (see 
jl]) if the probability that a random graph G n (p) has the property Q converges to 1 as n 
tends to infinity. Occasionally, we shall write almost all instead of almost every. It is easy 
to see that if F 1 is a pointmass at 0, i.e., Fi(x) = 1 for x > and F^x) = for x < 0, and 
i*2 is the Bernoulli distribution with mean p, then the Wigner matrix X n coincides with 
the adjacency matrix A(G n (p)) of the random graph G n (p). Obviously, o"2 = a/p(1 — p) 
in this case. 

To establish the exact estimate of the energy S'(G n (p)) for a.e. graph G n (p), we first 
present some notions. In what follows, we shall use A to denote the adjacency matrix 
A(G n (p)) for convenience. Set 

A = A-p(J n -I n ). 
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It is easy to check that each entry of A has mean 0. We define the energy <^(M) of a 
matrix M as the sum of absolute values of the eigenvalues of M. By virtue of the following 
two lemmas, we shall formulate an estimate of the energy S'(A), and then establish the 
exact estimate of S{A) = £{G n {p)) by using Lemma YTM 
Let I be the interval [—1,1]. 

Lemma 2.2. Let I c be the set R\I. Then 

lim / x 2 d^ n -i/ 2 -x(x) = / x 2 d&(x) a.s. 

n-^oo Jj c ■ Jj c 

Proof. Suppose 4> n -i/ 2 -^(x) is the density of $ n -i/a-x(x). According to Eq.flT]), with prob- 
ability 1, <f> n -i/3~\{%) converges to 4>(x) almost everywhere as n tends to infinity. Since 
4>{x) is bounded on /, it follows that with probability 1, x 2 4> n -i/2x( x ) is bounded almost 
everywhere on /. Invoking bounded convergence theorem yields 

lim / x 2 d^ n -i /2 -^(x) = / x 2 d&(x) a.s. 
n ^°° J i J i 

Combining the above fact with Eq.(j2J), we have 

lim / x 2 d$ n - l/2 - x (x) = lim f / 2 2 d$ n _i/ 2 x0) - / x 2 d$ n -i /2 - x (x) J 

n^oojj c n— oo \J Jj J 

= lim / x 2 d$ n -i/2-x(x) - lim / x 2 d^ n -i /2 -^(x) 

n— >oo J n— »oo J j 

= J x 2 d$(x) - J x 2 d§(x) a.s. 
x 2 d<&(x) a.s. 

□ 



Lemma 2.3 (Billingsley [3J pp. 219). Let fi be a measure. Suppose that functions 
a n , b n , f n converges almost everywhere to functions a, b, f , respectively, and that a n < f n < 
b n almost everywhere. If J a n d\i — > J adfi and J ' b n dfi — > jbd/u,, then J f n dfi — > J fdfi. 

We now turn to the estimate of the energy <?(A). To this end, we first investigate 
the convergence of J \x\d§ n -i /2 -^(x). According to Eq.flT]) and the bounded convergence 
theorem, we can deduce, by an argument similar to the first part of the proof of Lemma 
ET21 that 



lim / \x\d^ n -i/ 2 -^(x) = / |x|d$(a:) a.s. 

n~*oo J j J j 



Obviously, \x\ < x 2 if x G I c := R \ I. Set a n (x) = 0,b n (x) = x 2 (p n -i/ 2 -x(x), and 
f n (x) = \x\(p n -i/ 2 ^(x). Employing Lemmas 12.21 and ETBT we have 

lim / \x\d^ n -i /2 -^(x) = / a.s. 

n-*oo Jj c Jj c 

Consequently, 

lim / \x\d^ n -i/ 2 -^(x) = / |x|c/$(x) a.s. (3) 

n~*oo J J 



1 



Suppose Ai,...,A„ and A 1; ...,A n are the eigenvalues of A and n X / 2 A, respectively. 
Clearly, Y^i=i M = n1 ^ 2 J2i=i W\- By Eq.flSJ), we can deduce that 



<?(A)/" 3/2 = -^X>l 

i=l 

n 
i=l 

= J \x\d^ n - 1/is (x) 
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x\d<&(x) a.s. (n — > oo) 

2(72 



4 o"o —x 2 dx 



= (To = \/ p(l — p). 

3tt 3tt V ' PV ^ 
Therefore, the energy S (A) enjoys a.s. the equation as follows: 

* (A) = ( |- yW^P) + • 

We proceed to investigating S{A) = S{G n {p)) and present the following result due to 
Fan. 

Lemma 2.4 (Fan [6J). Let X, Y, Z &e rea/ symmetric matrices of order n such that 
X + Y = Z. Then 

n n n 

^Tia^i + ^imy^^Tia^z)! 

i=l i=l i=l 

where Aj(M) (i = 1, • • • , n) are the eigenvalues of the matrix M. 

It is not difficult to verify that the eigenvalues of the matrix J n — I n are n — 1 and 
— 1 of n — 1 times. Consequently <S{3 n — I n ) = 2(n — 1). One can readily see that 
S (p(J n -I re )) = p<f(J n - I re ). Thus, 

^(p(Jn-In)) =2p(n-l). 

Since A = A + p(J n — I n ), it follows from Lemma [2.41 that with probability 1, 
<f(A) < ^(A) + < f(p(J n -I n )) 



— ^(1-^ + 0(1)) +2p(n-l). 



Consequently, 



lim £(A)/n 3/2 < —y/ p (l-p) a.s. (4) 

n— >oo 37T 
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On the other hand, since A = A + p( — (J n — I n )) , we can deduce by Lemma 12.41 that 
with probability 1, 

£{A) > S (A) - (p( - (J n - I re ))) 

= S (A) - ^(p(J„ - I„)) 



— Vp(1-p) + o(1)J -2p(n-l). 

07T 



Consequently, 



lim <f(A)/n 3/2 > — ^/p(l - p) a.s. (5) 

n— >oo 37T 



Combining Ineq.(jl]) with Ineq.([5]), we have 



<g(A) = n^ [ y/rfT=d + (1) 

07T 



a.s. 



Recalling that A is the adjacency matrix of G n (p), we thus obtain that a.e. random graph 
G n (p) enjoys the equation as follows: 



Remark 2.2. Note that for p — |, Nikiforov in [12] got the above equation. Here, our 
result is for any probability p, which could be seen as a generalization of his result. 



3 The energy of the random multipartite graph 

We begin with the definition of the random multipartite graph. We use K n . Ult )Vm to 
denote the complete m-partite graph with vertex set [n] whose parts Vi, . . . , V m (m = 
m(n) > 2) are such that \Vi\ = nvi = ni/j(n), % — 1, . . . ,m. Let Q n ;v 1 ...v m ijP) be the set of 
random m-partite graphs with vertex set [n] in which the edges are chosen independently 
with probability p from the set of edges of K n . Vu ,^ Vm . We further introduce two classes 
of random m-partite graphs. Denote by G n ,m(p) an d G' nm (p), respectively, the sets of 
random m-partite graphs satisfy, respectively, the following conditions: 

lim max{ui(n) , . . . , v m (n)} > and lim % = 1. (6) 

and 

lim max{z/i(n), . . . , ^ m (n)} = 0. (7) 

n— >oo 

One can easily see that to obtain the estimate of the energy for the random multipartite 
graph G n]V1 _ Um {p) e Q n ;v\...u m {p)i we need to investigate the spectral distribution of the 
random matrix A(G n;vi _ Vm (p)). It is not difficult to verify that A(G n;vi _ Vm (p)) would 
be a special case of a random matrix X n (^i, . . . , z/ m ) (or X„^ m for short) called a random 
multipartite matrix which satisfies the following properties: 

independent random variables with 
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• the Xjj's have the same distribution Fi if i and j G while the Xi/s have the same 
distribution F 2 if z G and j G [n] \ Vjt, where Vi, . . . ,V m are the parts of i^n;^,...,^ 
and /c is an integer with 1 < < m; 

• l^ijl < K for some constant K. 

Apparently, if Fi is a pointmass at and F 2 is a Bernoulli distribution with mean p, then 
the random matrix X n m coincides with the adjacency matrix A(G n;J/1 ... l/m (p)). Thus, we 
can readily evaluate the energy S{G nMl ,,, Vm [p)) once we obtain the spectral distribution of 
X„ ]m . In fact, the random matrix X n m is a special case of the random matrix considered 
by Anderson and Zeitouni [1] in a rather general setting called the band matrix model 
which can be regarded as one of generalization of the Wigner matrix, and we shall employ 
their results to deal with the spectral distribution of X„ im . 

The rest of this section will be divided into two parts. In the first part, we shall present, 
respectively, exact estimates of the energies for random graphs G n>m (p) G G n ,m(p) an d 
G' n mip) e Q'n m(p) by exploring the spectral distribution of the band matrix. We establish 
lower and upper bounds of the energy for the random multipartite graph G n:ui __ Mm (p), 
and moreover we obtain an exact estimate of the energy for the random bipartite graph 
G n - UlyV2 {p) in the second part. 



3.1 The energy of G n>m (p) and G' nm {p) 

In this part, we shall formulate exact estimates of the energies for random graphs G n ^ m {p) 
and G' nm {p), respectively. For this purpose, we shall establish the following theorem. To 
state our result, we first present some notations. Let l n , m = (i Pt q) n xn be a quasi-unit 
matrix such that 



1 'dp, qeV k , 

if p G Vfc and q G [n] \ 14, 



where V\, . . . , V m are the parts of K n . Ult ^ Vm and k is an integer with 1 < k < m. Set 
fii = J xdFi (z = 1,2) and 

~^-n,m ~^-n,rn l^l^-n,m /^(J-n In,m)- 

Evidently, X n m is a random multipartite matrix as well in which each entry has mean 0. 
To make our statement concise, we define A 2 = (a 2 +(m — 1) cr|)/m. 

Theorem 3.1. 

(i) If condition fidj) holds, then 

i.e., the ESD $„-i/2x (x) converges weakly to a distribution ^>(x) in probability as 
n tends to infinity where ^{x) has the density 

^) = T^V^A 2 - X 2 1 W <2A. 
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(ii) If condition holds, then $ n -i/2x n m ( x ) as n — > oo. 

Our theorem can be proved by a result established by Anderson and Zeitouni pQ. We 
begin with a succinct introduction of the band matrix model defined by Anderson and 
Zeitouni in pQ, from which one can readily see that a random multipartite matrix is a 
band matrix. 

We fix a non-empty set C = {c\, C2, . . . , c m } which is finite or countably infinite. The 
elements of C are called colors. Let k be a surjection from [n] to the color set C, and 
we say that is the color of %. Naturally, we can obtain a partition Vi, . . . , V m of [n] 
according to the colors of its elements, i.e., two elements % and %' in [n] belong to the same 
part Vj if and only if their colors are identical. We next define the probability measure 
9 m on the color set as follows: 

e m (C) = e m{n) (C) = \K~\C)\/n, 1 < i < m = m(n), 

where C C C and k _1 (C) = {x G [n] : G C}. Evidently, the probability space 

(C,2 c ,9 m ) is a discrete probability space. Set 

e = iim e rn . 

For each positive integer we fix a bounded nonnegative function on color set 
and a symmetric bounded nonnegative function on the product of two copies of the 
color set. We make the following assumptions: 

• d^ is constant for k ^ 2; 

• s( k > is constant for k {2,4}; 

Let {Cij}i,j=i be a family of independent real- valued mean zero random variables. We 
suppose that for all 1 < i, j < n and positive integers k, 



E(ie«n < 



S ( fc )( K (z), K (j)) if^J, 

d^(K(i)) if i = J, 



and moreover we assume that equality holds above whenever one of the following condi- 
tions holds: 

• k = 2, 

• i j and k = 4. 

In other words, the rule is to enforce equality whenever the not- necessarily-constant func- 
tions d^ 2 \ or s^ 4 ) are involved, but otherwise merely to impose a bound. 

We are now ready to present the random symmetric matrix Y n called band matrix in 
which the entries are the r.v. Evidently, Y n is the same as X„ jm providing 

» <a) (*(0,«(j)) = {^ tfK(oJ«(j) ^d^\n{i)) = al^<^,J<n. (8) 



8 



So the random multipartite matrix X n m is a special case of the band matrix Y n . 

Define the standard semi-circle distribution $0,1 of zero mean and unit variance to be 
the measure on the real set of compact support with density 0o,i( x ) = — x 2 l| x |<2- 

Anderson and Zeitouni investigated the LSD of Y n and proved the following result (The- 
orem 3.5 in pQ). 

Lemma 3.2 (Anderson and Zeitouni // J s^ 2 \c,d)9{dd) = 1, then $ n -i/ 2Yn (x) 
converges weakly to the standard semi-circle distribution $ ,i ^ n probability as n tends to 
infinity. 

Remark 3.1. The main approach employed by Anderson and Zeitouni to prove the 
assertion is a combinatorial enumeration scheme for the different types of terms that con- 
tribute to the expectation of products of traces of powers of the matrices. It is worthwhile 
to point out that by an analogous method called moment approach one can readily ob- 
tain a stronger assertion for X n m that the convergence could be valid with probability 1. 
Moreover, one can show that for each positive integer k, 

x k ^/(x) a.s. if condition ([6]) holds, 
Jim I x k ^ 1/2Xn (x) = { J f (9) 

x $( s. if condition (T7D holds. 



However, we shall not present the proof of Eq.Q here since the arguments of the two 
methods are similar and the calculation of the moment approach is rather tedious. We 
refer the readers to Bai's survey [2] for details. 

Using Lemma [3.21 to prove Theorem 13.11 we just need to verify J s^ 2 \c, c')O(dc') = 1. 
For Theorem I3.1( i). we consider the matrix A _1 X n m where A 2 = (er 2 +(m — l)o~ 2 )/m. 
Note that condition ([6]) implies that 6 m {ci) ->l/m asn-^ 00, 1 < i < m. By means of 
condition (jSJ), one can readily see that for the random matrix A _1 X n m , 

A^ \m m J 

Consequently, Lemma [3.21 implies that 

$ n-i/2 A -ix n , m $o,i as n > 00. 

Therefore, 

$ n-i/2x„, m as n ^ 00, 

and thus the first part of Theorem 13.11 follows. 

For the second part of Theorem 13. 1\ we consider the matrix 1 X n m . Note that 
condition ((7|) implies that 9(cj) = lim^oo 9 m (ci) = lim^oo z/j(n) = 0, 1 < i < m. By 
virtue of condition (jBJ), if c 7^ d then s^ 2 \c, d) = 1. Consequently, for the random matrix 
ao 1 X„ m , we have 

W(c,d)9(dd) = [ sW(c,c') Xc , {c} 9(dd) 



X cx{c) 0{dd) 
»(C\{c}) = l. 
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As a result, Lemma [3.21 implies that 

$„-i/2.-iv $0,1 as n — > oo. 

Therefore, 

$ „-V2x„, m $0*0 as n -> oo, 

and thus the second part follows. 

We now employ Theorem 13. II to estimate the energy of G n . Ul ^ Mm (p) under condition 
or (J7|). For convenience, we shall use A nim to denote the adjacency matrix A.{G n>m {p)). 
One can readily see that if a random multipartite matrix X n m satisfies condition ©, and 
Fx is a pointmass at and F2 is a Bernoulli distribution with mean p, then X„ jm coincides 
with the adjacency matrix A nj?n . Set 

A}2 jm A nm p(3 n In,m) (^-^) 

where I„ ;m is the quasi-unit matrix whose parts are the same as A n)Jn . Evidently, each 
entry of A n m has mean 0. It follows from the first part of Theorem 13. II that 

$.-i/2T ^(x) as n — > 00. 

Since the density of ty(x) is bounded with the finite support, we can use a similar method 
for showing Eq.(j3J) to prove that 

j \ x \ d ®n-v*A n , m ( x ) J \x\d&(x) as n -> 00. 



Consequently, 

^(A n , m )/n 3 / 2 = J \x\d$ n - 1/2Xn Jx) 

— >p J \x\d^(x) as n — > 00 



m /"V^?^ / (m-l)o-i ~ , 
xh/4 or dx 



27r(m — 1) (Tn J _2 x pnEl 

V m 



m — 1 8 /m — 1 

1/ °"2 = P(l - P • 

37T V 571 37T V 571 

Therefore, a.e. random matrix A njTTl enjoys the equation as follows: 



* (A,, m ) = n 3 / 2 f^V P(1_p) + ° (i; 
We now turn to the estimate of the energy #(A rejm ) = S{G n>m {p)). Evidently, 

Jn In,m (Jn In) "I - (In In,m)- 

By virtue of Lemma I2.4[ we have 

^(Jn — In,m) — ^(Jn — In) + ^(In — In,m)- 
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Recalling the definition of the quasi-unit matrix I n m and the fact that <S{3 n — I n ) = 
2(n — 1), we have <f(J„ — I n ,m) < 0(n). According to Eq.l TTUl) , we can use a similar 
argument for the estimate of the energy <?(A) from ^(A) to show that a.e. random 
matrix A n m enjoys the equation as follows: 



* (A n , m ) = ± J —p(l - P ) + o(l) . 

\ 3n V m I 

Since the random matrix A n , jm is the adjacency matrix of G n ,m(p), we thus show that a.e. 
random graph G n>m {p) enjoys the following equation: 



In what follows, we shall use A! n>m to denote the adjacency matrix A(G' nm (p)). It is 
easily seen that if a random multipartite matrix X nm satisfies condition ([71) , and Fi is a 
pointmass at and F 2 is a Bernoulli distribution with mean p, then X n m coincides with 
the adjacency matrix A' n m . Set 

A' ' n ,m ■^-n,m Pi^n ^-n,m) 

where I' nm is the quasi-unit matrix whose parts are the same as A' nm . One can readily 
check that each entry in A' n m has mean 0. It follows from the second part of Theorem 
EI that 

$ B -W^( Z ) as n oo. 

Employing the argument analogous to the estimate of $(p{3 n — I n ,m)), <^(A„ im ) and 
^(A„ jm ), one can evaluate, respectively, $(p{3 n — I' nm )), <f(A' njm ) and $(A! nm ), and 
finally show that a.e. random graph G' nm (p) satisfying condition (J7J enjoys the following 
equation: 

*(Gn,M) = n 3/2 (J^VpO^P) + o(l)) • (11) 
3.2 The energy of G n]Vl .. Mm (p) 

In this part, we shall give an estimate of energy for the random multipartite graph 
G n;Ul .. Mm (p) satisfying the following condition: 

lim max{z/i(n), . . . , v m {n)} > and there exist v\ and Uj , lim 1 < 1. (12) 

Moreover, for random bipartite graphs G n;vi)V2 {p) satisfying \im n ^ ^ Ui(n) > (i = 1,2), 
we shall formulate an exact estimate of its energy. 

Anderson and Zeitouni [1] established the existence of the LSD of X n>m with partitions 
satisfying condition (|12[) . Unfortunately, they failed to get the exact form of the LSD, 
which appears to be much hard and complicated. However, we can establish the lower 
and upper bounds for the energy <S{G n . Vl ^ Mm {p)) via another way. 
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Here, we still denote the adjacency matrix of multipartite graph satisfying condition 
( fl2l) by A.„ jm . Without loss of generality, we assume, for some r > 1, |Vi|, . . . , \ V r \ are of 
order 0{n) while |K+i|, • ■ • , \V m \ of order o{n). Let A' nm be a random symmetric matrix 
such that 

(ij) if % or j £ V s , 1 < s < r, 
tij if i,j e V s , 1 < s < r and i > j, 

if i, j E V s (r + 1 < s < m) or i = j, 

where ti/s are independent Bernulli r.v. with mean p. Evidently, A' nm is a random mul- 
tipartite matrix. By means of Eq. lflTT) . we have S{A' nrn ) = (^^p(l — p) + o(l)j n 3//2 . 



Set 



r> — a' — a 

" n m, Jr ^-n,m 



\ 







(13) 



V 0/ 

One can readily see that K,; (z = 1, . . . , r) is a Wigner matrix and thus a.e. Kj enjoys the 
following: 



*(K,) = f -v^T^) + (l) ) {v in f 2 . 
Consequently, a.e. matrix D n satisfies the following: 



<f(D n ) = ( — vMl-J>) + o(l) ) (v 1 2 +--- + z/ r 5 )n 



By Eq.([I3]), we have A n>m + D n = A' nm and A' n m + (-D n ) = A n>m . Employing Lemma 
12.41 we deduce 

£(A! n ,m) - *(D B ) < ^(A n , m ) < *(A^ J + *(D B ). 
Recalling that A n m is the adjacency matrix of G n . Ul ^ Mm (p), the following result is relevant. 
Theorem 3.3. Almost every random graph G n . ul .__ Um (p) satisfies the inequality below 



^ i,^< < f(G„ 11 , I ...„,„(p))^v5a :r rt+o(i)j < (i + E 



-i 



3/2 



Remark 3.2. Since ui,...,u r are positive real numbers with 2^[=i v i — 1> we have 

YZ=i v i(l ~ ^ 1/2 ) > °- Therefore, Yh=\ u % > Yfi=i u i /2 ^ and tnus 1 > SLi^f ■ Hence, 
we can deduce, by the theorem above, that a.e. random graph G nw .„ % (p) enjoys the 
following 

<?(G n ^ m (p)) = 0(nV 2 ). 
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In what follows, we investigate the energy of random bipartite graphs G n]UltU2 (p) sat- 
isfying lim^oo Vi{n) > (i — 1,2), and present the precise estimate of &{G n - Vl}Vi (p)) via 
Marcenko-Pastur Law. 

For convenience, set ri\ = v\n and n 2 = v 2 n. Let be a quasi-unit matrix with the 
same partition as A n2 - Set 



^n,2 — A ni 2 — p(3 n — In,2) 



o X T 
X o 



(14) 



where X is a random matrix of order n 2 x n\ in which the entries X(ij) are iid. with 
mean zero and variance a 2 = p(l — p). By the equation 

AI ni \f AI ni -X T \ _ \ f -X T 

-X AI„ 2 J V AI,„ 2 - \- 1 XX T J V -X AI n2 



we have 

A n ■ A ni |AI n2 — A 1 XX r | = A™|AI n — A nj2 | 



and consequently, 

A ni |A 2 I„ 2 — XX T | = A" 2 |AI„ — A„ j2 | 



Thus, the eigenvalues of A nj 2 are symmetric, and moreover A is the eigenvalue of -j= A nj 2 

if and only if A is the eigenvalue of ^-XX T . Therefore, we can characterize the spectral 
of A n .2 by the spectral of XX T . Bai formulated the LSD of ^-XX T (Theorem 2.5 in [2]) 
by moment approach. 

Lemma 3.4 (Marcenko-Pastur Law [2]). Suppose that X(zj) 's are iid. with mean 
zero and variance a 2 = p(l —p), and v 2 jv\ — * V £ (0, oo). Then, with probability 1, the 
ESD $i_ XX T converges weakly to the Marcenko-Pastur Law F y as n — > oo where F y has 
the density 

f y (x) = ^.^y Vib-^x-a) l a <,< fe 

and has a point mass 1 — 1/y at the origin if y > 1 where a = p(l — p)(l — y^) 2 and 
b = p(l-p)(l + jy) 2 . 

By the symmetry of the eigenvalues of -^==A nj2 , to evaluate the energy $ (-^=A re)2 ), 



J2 1- 

we just need to consider the positive eigenvalues. Define Q n2 (x) = — ^ 2L - One can 
see that the sum of the positive eigenvalues of -4=A nj2 equals n 2 J °° xdO n2 (x). Suppose 
< X\ < x 2 , we have 

Qn 2 (x 2 ) - e n2 {xx) = $±xXt(^) - ^XX^l)- 

It follows that 



roo roc 

/ xdB n2 (x) = / \fxd<& 
Jo Jo 



13 



Note that all the eigenvalues of ^-XX T are nonnegative. By the moment approach (see 
[2] for instance), we have 

/poo 
X 2 d$^ xxT (x) = / X 2 d$^ xxT (x) 

POO 

— > / x 2 dF y (x) a.s. (n — > oo) 
x 2 c?Fj / (x) 

Analogous to the proof of Eq.([3]), we can deduce that 

POO POO 

lim / v / xrf ( l ) j_ XX T(x) = / \fxdF y (x) a.s. 

n^ooJ Q n-i ■ J 



Therefore, 

"oo rVb 

7rp(l -p)y 



POO PVU 

lim / xdQ n2 {x) = / — — <J (b — x 2 ){x 2 — a) dx a.s. 



Let 

A= / — -a/^-x 2 )^ 2 -^ dx. 

iv^ np(l-p)y 

We obtain that for a.e. A nj2 the sum of the positive eigenvalues is (A + o(l))n2y / nT. Thus, 
a.e. <^(A raj2 ) enjoys the equation as follows: 



^(A„ )2 ) = (2A + o(l))n 2v /^. 

Furthermore, we can get 

Vb[(a + 6) Ep(l - a/b) - 2a Ek(l - a/b)} 



A 



37rp(l 



where Ek is the complete elliptic integral of the first kind and Ep is the complete elliptic 
integral of the second kind. Let t £ [0, 1], the two kinds of complete elliptic integral are 
defined as follows 

r% d9 f% i 

Ek(t) = / - and Ep(t) = / \/l-tsm 2 6 d6. 

Jo y/l-t sin 2 9 Jo 

The value can be got by mathematical software for every parameter t. 
Employing Eq. (ll4p and Lemma [2.41 we have 

S (A„, 2 ) - <?(p(J n - l n>2 )) < #(A ni2 ) < <f (A n , 2 ) + £(j>{J n - I n , 2 )). 

Together with the fact that <o(p(J n — l n ,2)) = 2py / Z^~z/ 2 ~n and n 2 ^/n\ = ^ 2 \/^i^ 3 ^ 2 , we get 

# (A n , 2 ) = (2^ 2V /^A + o(l))n 3 / 2 . 

Therefore, the following theorem is relevant. 
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Theorem 3.5. Almost every random bipartite graph G n]UltU2 (p) with z^/^i ~~ ► V enjoys 

&(GwM) = (2^ 2v ^A + o(l))r i 3 / 2 . 

We now compare the above estimate of the energy ${G n . vl>U2 {p)) with bounds obtained 
in Theorem 13.31 for p = 1/2. In fact, Koolen and Moulton [11] established the following 
upper bound of the energy ${G) for simple graphs G: 

71 

^(G)<-(V^+1). 

Consequently, for p = 1/2, this upper bound is better than ours. So we turn our attention 
to compare the estimate of S{G n -^ XyV2 (l/2)) in Theorem 13 . 51 with the lower bound in Theo- 
rem 13.31 By the numerical computation (see the table below), the energy £'(G n . Ul!U2 (l/2)) 
of a.e. random bipartite G n;vijl/2 (l/2) is close to our lower bound. 



y 


S{G n . vl , V2 {p)) 


lower bound of $(G n ^ V2 {p)) 


l 


(0.3001 + o(l))n 3/2 


(0.1243 + o(l))n 3 / 2 


2 


(0.2539 + o(l))n 3 / 2 


(0.1118 + o(l))n 3 / 2 


3 


(0.2071 + o(l))n 3 / 2 


(0.0957 + o(l))n 3 / 2 


4 


(0.1731 + o(l))n 3 / 2 


(0.0828 + o(l))n 3 / 2 


5 


(0.1482 + o(l))n 3 / 2 


(0.0727 + o(l))n 3 / 2 


6 


(0.1294 + o(l))n 3 / 2 


(0.06470 + o(l))n 3 / 2 


7 


(0.1148 + o(l))n 3 / 2 


(0.05828 + o(l))n 3 / 2 


8 


(0.1031 + o(l))n 3 / 2 


(0.05301 + o(l))n 3 / 2 


9 


(0.09353 + o(l))n 3 / 2 


(0.04862 + o(l))n 3 / 2 


10 


(0.08558 + o(l))n 3 / 2 


(0.04491 + o(l))n 3 / 2 
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